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Abstract Supermassive black holes (SMBHs) are common in local galactic nuclei, 

and SMBHs as massive as several billion solar masses already exist at redshift z = 6. 

Q_^ . These earliest SMBHs may grow by the combination of radiation-pressure-limited 

Q^ ' accretion and mergers of stellar-mass seed BHs, left behind by the first generation 

P-n . of metal-free stars, or may be formed by more rapid direct collapse of gas in rare 

special environments where dense gas can accumulate without first fragmenting into 
stars. This chapter offers a review of these two competing scenarios, as well as some 
more exotic alternative ideas. It also briefly discusses how the different models may 
£i , be distinguished in the future by observations with JWST, LISA and other instru- 

ct ' ments. 
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1 Introduction 



The discovery of about two dozen bright quasars with luminosities > 10 47 erg s _1 
at redshift z — 6 suggests that some supermassive black (SMBHs) as massive as a 
fewx 10 9 M Q have been already assembled when the universe was less than 1 Gyr 
old (see, e.g., ref. AMI for a review). These high-redshift quasars are exceedingly 
rare, with a space density of order ~ lGpc~ 3 , and can only be found in large sur- 
veys of the sky, such as the Sloan Digital Sky Survey (SDSS), or the smaller-area 
but deeper CFHQS [284] and UKIDSS 1 14611 surveys. These quasars overall appear 
to be "fully developed", with spectra and metallicity patterns that appear remark- 
*-j ■ ably similar to their counterparts at moderate redshifts [|62l . Indeed, if one selects 

individual quasars with the same luminosity, their properties show little evolution 
with cosmic epochjj This implies that the behavior of individual quasars is proba- 
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bly determined by local physics near the SMBH and is not directly coupled to the 
cosmological context in which the SMBH is embedded. However, it is clear that 
the quasar population as a whole does evolve over cosmic timescales. Observations 
from < z < 6 in the optical (e.g., the Anglo-Australian Telescope's Two Degree 
Field, or 2dF, and the Sloan Digital Sky Survey, or SDSS) and radio bands 12281 
show a pronounced peak in the abundance of bright quasars at z ~ 2.5. A similar 
behavior has been confirmed in X-ray observations [12381 . 

The cosmic evolution of quasar black holes between < z < 6 is likely driven 
by a mechanism other than local physics near the hole. This is reinforced by the 
fact that the timescale of activity of individual quasars is significantly shorter than 
cosmic timescales at z < 6, both on theoretical grounds (~4x 10 7 yr, the e-folding 
time for the growth of mass in a SMBH, whose accretion converts mass to radiation 
with an efficiency of e = Mc 2 /LEdd ~ 10%) and is limited by its own [Eddington] 
luminosity), and using the duty cycle of quasar activity inferred from various obser- 
vations (also ~ 10 7 yr but with large uncertainties, e.g. 01601 and references therein; 
see also |[9"T1 l2l8ll2"ITll2"22l ). 

In the cosmological context, it is tempting to link the evolution of massive quasar 
black holes with that of dark matter (DM) halos condensing in a Cold Dark Mat- 
ter (CDM) dominated universe, as the halo population naturally evolves on cosmic 
timescales l57l . Indeed, this connection has proven enormously fruitful and has re- 
sulted in the following broad picture: the first massive astrophysical black holes ap- 
pear at high redshifts (z > 10) in the shallow potential wells of low mass (< 10 8 M©) 
dark matter halos. These black holes grow by mergers and gas accretion, evolve into 
the population of bright quasars observed at lower redshifts, and eventually leave the 
SMBH remnants that are ubiquitous at the centers of galaxies in the nearby universe. 

In this picture, the presence of fewx 10 9 M© SMBHs at z > 6 presents a puz- 
zle j99o Metal-free stars, with masses ~ 100 M Q , are expected to form at red- 
shifts as high as z > 25 131 1351 12991 . and leave behind remnant BHs with similar 
masses f38ljl06| . However, the natural time-scale, i.e. the Eddington time, for grow- 
ing these seed BHs by > 7 orders of magnitude in mass is comparable to the age 
of the universe (e.g. ref. 11991 ). This makes it difficult to reach 10 9 M© without a 
phase of rapid (at least modestly super-Eddington) accretion, unless a list of opti- 
mistic assumptions are made in hierarchical merger models, in which multiple seed 
BHs are allowed to grow without interruption, and to combine into a single SMBH 
Il86l 12961 [331 12261 12761 Il9llll50l 12351 I25U1 . 

An alternative class of explanations involves yet more rapid gas accretion or col- 
lapse HHl [37] d38] Q32] [240] |25] EH] [2871 ESS EE1 EH) ■ In this family of models, 
primordial gas collapses rapidly into a SMBH as massive as 10 4 — 10 6 M©, either 
directly, or possibly by accreting onto a pre-existing smaller seed BH 112751 . or go- 
ing through the intermediate state of a very massive star 11371 . a "quasistar" l23l . or 
a dense stellar cluster 1186] [49] ■ These so-called "direct collapse" models involve 
metal-free gas in relatively massive (> 10 8 M©) dark matter halos at redshift z > 10, 



2 More generally, the non-trivial cosmological implications of the existence of massive BHs at 
early times was noted already when quasars were first found at redshifts z > 4 12591 . 
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with virial temperatures r v ; r > 10 4 K. The gas that cools and collapses in these halos 
must avoid fragmentation, shed angular momentum efficiently, and collapse rapidly. 

Many uncertainties about each of the above scenarios remain, and the astrophysi- 
cal process(es) responsible for the formation of the earliest massive black holes (and 
indeed for the presence of SMBHs at all redshifts) remain poorly understood. In 
this review, we focus on the emergence of the first generation of black holes, though 
many of the important questions are quite general and apply equally to subsequent 
generations of black holes. This review is organized as follows. In § [23 we describe 
theoretical expectations for the formation and growth of these black holes within the 
paradigm of hierarchical CDM cosmologies. In §|3] we "zoom in" and consider the 
local physics of black hole formation, and various pathways which could lead to the 
early presence of supermassive black holes. In §|4] we summarize several relevant 
recent observations that have implications for early black holes, and speculate on 
the power of future observations to probe the physics of the first BHs. We offer our 
conclusions in §|5] 

This chapter is an expanded an updated version of an earlier review 01001 . An- 
other recent, complimentary review on SMBH formation at high redshift can be 
found in J27T]. 



2 First Structure Formation 

In this section, we sketch some basic theoretical arguments relevant to the forma- 
tion of structure in the early universe. We then discuss formation mechanisms for 
SMBHs. 



2.1 Cosmological Perturbations as the Sites of the First Black 
Holes 

Measurements of the Cosmic Microwave Background (CMB) anisotropies by the 
Wilkinson Microwave Anisotropy Probe (WMAP), determinations of the luminosity 
distance to distant type la Supernovae, and other observations have led to the emer- 
gence of a robust "best fit" cosmological model with energy densities in CDM and 
"dark energy" of (Q M ,Q A ) « (0.3,0.7) (see, e.g. Ifl37ll , for the seven-year WMAP 
results, and its combination with other datasets). 

The growth of density fluctuations and their evolution into nonlinear dark mat- 
ter structures can be followed in this cosmological model from first principles by 
semi-analytic methods II 193112301 . More recently, it has become possible to derive 
accurate dark matter halo mass functions directly in large cosmological N-body sim- 
ulations HI 181 . with different codes agreeing at the 10% level, and mass functions 
measured down masses as low as ~ 10 6 Mq and high redshifts as high as z ~ 30 
(e.g. 1153J |2001). 
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Within the ACDM model, with a scale-invariant primordial power spectrum, ro- 
bust predictions can therefore be made for the dark matter halos. Structure formation 
in such a universe is "bottom-up", with low-mass halos condensing first. Halos with 
the masses of globular clusters (10 5 ~ 6 M ) are predicted to have condensed from 
~ 3<T peaks of the initial primordial density field as early as ~ 1% of the current age 
of the universe, or at redshifts of z ~ 25. These predictions are limited mainly by 
the 5 — 10% uncertainty in the normalization of the primordial power spectrum, as, 
and by the need to extrapolate the power-spectrum 2-3 orders of magnitude below 
the scales on which it has been directly constrained. In warm dark matter models, 
with particle masses of order ~ IkeV or less, free-streaming would result in a sharp 
exponential suppression of the fluctuation power on the relevant scales (masses be- 
low 10 8 M ), and could significantly reduce the number of DM halos at the earliest 
redshifts |fl3l 13001 . 

It is natural to identify the first collapsed DM halos as the sites where the first 
astrophysical objects, including the first black holes, were born. The nature of the 
objects that form in these early dark matter halos is currently one of the most rapidly 
evolving research frontiers in astronomy. 



2.2 Chemistry and Gas Cooling at High Redshifts 

Baryonic gas that falls into the earliest nonlinear dark matter halos is unable to cool 
efficiently, and is shock heated to the characteristic virial temperatures less than a 
few hundred Kelvin. It has long been pointed out [26, 202, 282| that such gas needs 
to lose its thermal energy efficiently (within about a dynamical time) in order to 
continue contracting, or in order to fragment. In the absence of any dissipation, it 
would simply reach hydrostatic equilibrium and would eventually be incorporated 
into a more massive halo further down the halo merger hierarchy. While the forma- 
tion of nonlinear dark matter halos can be followed from first principles, the cooling 
and contraction of the baryons, and the ultimate formation of stars or black holes in 
these halos, is much more difficult to model ab initio. 

The gas content of a cosmological perturbation can contract together with the 
dark matter only in dark halos above the cosmological Jeans mass, Mj w 10 4 M [(1 + 
z)/l l] 3 ' 2 , in which the gravity of dark matter can overwhelm thermal gas pressure. 
Recent work 112561 has shown that immediately following recombination (at redshift 
Z ~ 1 , 000), the baryons develop coherent streaming motions relative to the dark mat- 
ter, with relative speeds of order 30 km s , on scales of a few Mpc. These relative 
velocities decay as « (1 +z), and reduce to ~ 1 km s _1 by z ~ 30, comparable to 
the velocity dispersions in the smallest dark matter halos at this epoch. In the some- 
what more massive halos in which the baryons can typically cool efficiently (with 
velocity dispersions of several km s~ ' at z < 30; see below), the streaming motions 
are expected to have at most a modest effect, reducing the gas fraction within the 
virial radius by a factor of ~4wo [2551 1 1 59l 12441 [781 |48l . However, as we will see 
below, the "stellar-seed" model for SMBH growth relies on one (or a few) "special" 



The Formation of the First Massive Black Holes 5 

rare seed BHs that form in few x ~ 10 5 M Q halos at redshifts as high as z ~ 30. 
This scenario appears vulnerable to the streaming motions. 

In the earliest, chemically pristine clouds, radiative cooling is dominated by H2 
molecules. As a result, gas phase H2 "astro-chemistry" is likely to determine the 
epoch when the first stars and black holes appear (primordial molecular chemistry, 
focusing on the role of H2 early structure formation was reviewed by []4]). Several 
papers have constructed complete gas-phase reaction networks and identified the 
two possible ways of gas-phase formation of H2 via the H or H^ channels. These 
were applied to derive the H2 abundance under densities and temperatures expected 
in collapsing high redshift objects HUH] [J62l [190J HH] |224j [J26l [J25] [2231 . Studies 
that incorporate H2 chemistry into cosmological models and that address issues such 
as non-equilibrium chemistry, dynamics, or radiative transfer have appeared rela- 
tively more recently. Ref 01031 used spherically symmetric simulations to study the 
masses and redshifts of the earliest objects that can collapse and cool via H2; their 
findings were confirmed by a semi-analytic treatment |253]. The first fully three di- 
mensional cosmological simulations that incorporate H2 chemistry and cooling date 
back to refs HHED and Q. 

The basic picture that emerged from these papers is as follows. The H2 fraction 
after recombination in the smooth "protogalactic" gas is small (jch2 == "H2/"H ~ 
10~ 6 ). At high redshifts (z > 100), H2 formation is inhibited, even in overdense 
regions, because the required intermediaries Hj and H are dissociated by cos- 
mic "microwave" background (CMB, but with the typical wavelength then in the 
infrared) photons. However, at lower redshifts, when the CMB photons redshift to 
lower energies, the intermediaries survive, and a sufficiently large H2 abundance 
builds up inside collapsed clouds (xh2 ~ 10~ 3 ) at redshifts z < 100 to cause cooling 
on a timescale shorter than the dynamical time. Sufficient H2 formation and cool- 
ing is, however, possible only if the gas reaches temperatures in excess of ~ 200 K 
or masses of a few x 10 5 M©[(1 + z)/ll]~ 3 ' 2 (note that while the cosmological 
Jeans mass increases with redshift, the mass corresponding to the cooling threshold, 
which is well approximated by a fixed virial temperature, has the opposite behavior 
and decreases at high redshift). The efficient gas cooling in these halos suggests that 
the first nonlinear objects in the universe were born inside ~ 10 5 M dark matter 
halos at redshifts of z ~ 20 — 30, corresponding to an ~ 3 — 4(7 peaks of the primor- 
dial density peak (of course, yet rarer low-mass halos exist even earlier - the first 
one within our Hubble volume collapsing as early as z ~ 60 117611 ). 

The behavior of metal-free gas in such a cosmological "minihalo" is a well posed 
problem that has been addressed in three dimensional numerical simulations. The 
first series of such simulations ll2l l3l [34] l35l 12981 were able to follow the contraction 
of gas to much higher densities than preceding studies. They have shown conver- 
gence toward a temperature/density regime of T ~ 200 K, n ~ 10 4 cm~ 3 , dictated 
by the critical density at which the excited states of H2 reach equilibrium and cool- 
ing becomes less efficient ||67l . These simulations suggested that the gas does not 
fragment further into clumps below sizes of 10 2 — 10 3 M Q , but rather it forms unusu- 
ally massive stars. Very recent simulations reached higher resolution than the earlier 
ones, and, in some cases, using sink particles, were able to continue their runs be- 
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yond the point at which the first ultra-dense clump develops [258 2461 1771 [1951 . 
These simulations suggest that at least in some cases, the gas in the central regions 
does, eventually, fragment into two or more distinct clumps, raising the possibility 
that the first stars formed in pairs, or even in higher-multiple systems. 

The masses of the first stars would then presumably be reduced. The initial mass 
function (IMF) of the first stars is crucial, and is indeed one of the most impor- 
tant uncertainties for early BH formation. This is because massive stars would 
naturally leave behind black hole seeds, which can subsequently grow by merg- 
ers and accretion into the SMBHs. Interestingly, massive stars appear to have an 
"either/or" behavior. Non-rotating stars with masses between ~ 40 — 140 M© and 
above ~ 260 M fl collapse directly into a black hole without an explosion, and hence 
without ejecting their metal yields into the surrounding medium, whereas stars in 
the range ~ 140 — 260 M Q explode as pair-instability supernovae without leav- 
ing a remnant 01061 . In contrast, stars with initial masses ~ 25 — 40 M still leave 
BH remnants but also eject metals, whereas those with masses M < 25 M Q do not 
leave any BH remnants. This dichotomy is especially interesting because early mas- 
sive stars are attractive candidates for polluting the IGM with metals at high red- 
shifts [ 155, 28 1 j. It is likely that the first stars had a range of masses, in which case 
they could contribute to both metal enrichment and to the seed black hole popula- 
tion, with a relative fraction that depends sensitively on their initial mass function 
(IMF). 



3 Massive Black Hole Formation 

Having reviewed the general problem of structure formation at high redshifts, we 
now focus on the question of how the first SMBHs were assembled. It is worth 
emphasizing that this is an unsolved problem - indeed, it is not entirely clear even 
whether the first nonlinear objects in the universe were stars or black holes, and 
whether galaxies or their central black holes formed first [1871 . The leading ideas re- 
lated to the formation of SMBHs at high redshifts can be broadly divided into three 
areas: (1) formation of seed black holes from "normal" stellar evolution and sub- 
sequent Eddington-limited accretion, (2) rapid direct collapse of gas to a SMBH, 
usually via a supermassive star/disk, and (3) formation of a SMBH (or an IMBH 
seed) by stellar dynamical processes in dense stellar systems, such as star clusters 
or galactic nuclei. It is, of course, possible that all of these processes could be rele- 
vant 1221I20T1 . 
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3.1 Growth from Stellar-Mass Seeds 

3.1.1 Basic Ingredients and Uncertainties 

Perhaps the most natural possibility is that early SMBHs grow by gas accretion out 
of stellar-mass seed black holes, left behind by early generation of massive stars. If 
the subsequent gas accretion obeys the Eddington limit and the hole shines with a 
radiative efficiency of 10%, then the time it takes for a SMBH to grow to the size 
of 3 x 10 9 M from a stellar seed of - 100 M is 3 x 10 7 ln(3 x 10 9 /l00) yr ~ 
7 x 10 8 yr. This is comparable to the age of the universe at z = 6 (~ 9 x 10 8 yr 
for a flat ACDM universe with H = 70 km s _1 Mpc -1 and Q M = 0.3). Therefore, 
the presence of these SMBHs is consistent with the simplest model for black hole 
growth, provided that (i) the seeds are present early on (at z > 15; see below), (ii) 
and the near-Eddington growth is uninterrupted. As the ~ 10 M Q host halo of 
the initial seed BH gets incorporated into the ~ 10 12 ~ 13 M© host halo of the z ~ 6 
SMBH, it grows by ~ 7 orders in magnitude, and experiences a large number of 
mergers with other, comparable-sized, halos. These merger partners may (or may 
not) have a growing BH at their centers. Therefore, these mergers need to be taken 
into account, and the "stellar seed" model most likely can not be viewed as that of a 
single seed BH, growing in isolation. 

Several authors have worked out the growth of SMBHs from stellar-mass seeds, 
by following the build-up of dark DM halos, and using simple prescriptions to track 
the formation of seed BHs, their subsequent growth by accretion, and their mergers. 
This can be done either semi-analytically [99 295 86 226], using Monte-Carlo 
realizations of the DM merger trees H2961 1331 12761 12501 . or based on cosmologi- 
cal hydrodynamics simulations 01501 11911 12351 . As noted in the Introduction, the 
uncertainties about the statistics of the DM halo merger trees are essentially negligi- 
ble since DM halo formation has been directly resolved in numerical simulations 
at the relevant low masses (down to ~ 10 6 M©) and high redshifts (out to z « 30). 
The most important - and still highly uncertain - ingredients of this 'stellar seed' 
scenario can be summarized as follows. 

(i) What is the smallest possible mass (or virial temperature, T see d)for early DM 
halos in which PopIII stars can form? A reasonable answer is r see d =few x 100 K, 
which allows molecular H2-cooling 1103112531 . 

( ii) In what fraction (/ see dJ of these halos do seed BHs actually form? This is a 
much more difficult question, since various feedback processes (due to radiation, 
metal pollution, or mechanical energy deposition) could suppress PopIII star for- 
mation in the vast majority of early low-mass halos ( |88, 92 1; see also a recent 
review 11361 ). Interestingly, the WMAP measurement of the electron scattering opti- 
cal depth provides empirical evidence that such negative feedback took place early 
on and shaped the reionization history J89l . The answer also depends on the IMF of 



3 At least in principle, since halo mass functions in large N-body simulations agree at the few 
percent level. In practice, however, there can be significant disagreements between Monte-Carlo 



halo merger trees made with different algorithms 113051 . 
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PopIII stars, since, as noted above, whether the stars leave a BH remnant or explode 
as pair instability SNe depends on their masses. The dividing mass, rj 25 M , was 
evaluated in non-rotating stellar evolution models 01061 , whereas recent simulations 
indicate that the first stars in minihalos have significant rotation 12451 . Rotation can 
help drive winds and prevent BH formation entirely, or can produce a hypernova and 
reduce the mass of the remnant BH. Finally the velocity dispersions of the lowest- 
mass minihalos are only a few km s - 1 , only a factor of few higher than the residual 
bulk streaming motions between the gas and the DM halo at z > 20. These stream- 
ing motions can therefore reduce the gas fractions in the earliest minihalos and also 
lower the stellar masses by driving turbulence ll78l . 

(Hi) What is the time-averaged accretion rate of the seed BHs? This is con- 
veniently parameterized by a duty cycle /duty, defined as the fraction of the mass 
accretion rate that would produce the Eddington luminosity, if e ps 10% of the rest 
mass was converted to radiation (so that /d uty = 1 is the fiducial Eddington rate). 
Radiative feedback is usually expected to lead to sub-Eddington rates (e.g. J6|), 
and in spherical symmetry, the accretion was recently to shown to be episodic, with 
/duty ~ 0.3 11731 . The expectation is therefore that /j u ty is less than unity. In practice, 
if the accretion is radiatively inefficient, or if the radiation is trapped or is beamed 
and "leaks out", then /duty could exceed unity (see more on this below). 

(iv) Finally, what happens when DM halos merge? The simplest and most op- 
timistic assumption is that the BHs promptly coalesce, as well. However, even if 
dynamical friction on the DM (and on any stars present in later stages of the merger 
hierarchy) is efficient, it is possible that, due to the radiation of its parent star, the 
remnant BHs are no longer embedded in dense enough gas to allow this. Further- 
more, even if the BHs coalesce, the merged binary BH can suffer strong gravita- 
tional recoil at the time of the merger, due to the linear momentum carried by the 
anisotropic emission of gravitational waves (e.g. fTT) and references therein). Such 
a "kick" can eject the BH from the shallow potential wells (~ lkm/s) of the early 
halos, and the BH will be effectively lost. While kicks for comparable-mass BHs 
with random spins are of the order of ~ 100 km/s, the kick speed depends strongly 
on the mass ratio and on the spin vectors of the two BHs. In particular, kicks be- 
come very small (< lkm/s) for mass ratios q = M\ /M2 < 10~ 2 , irrespective of BH 
spins (e.g. ifTTI '). This may be key to avoid loosing growing seed BHs by ejection, 
and thus for the buildup of SMBHs early on. 



3.1.2 Worked Illustrative Examples 

In FigureQ] we show SMBH mass functions at z = 6, illustrating the impact of three 
of the most uncertain basic assumptions above, taken from a recent example of 
the Monte Carlo merger tree approach [250 1. The mass functions were constructed 
from the merger histories of w 10 5 DM halos with masses M > 1O 8 M0 at redshift 
z = 6. The upper right region in each panel, demarcated by the red dashed rectangle, 
shows the observational constraint on the SMBH space density. Each galaxy was 
modeled with a spherically symmetric mass distribution consisting of a DM halo 
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Fig. 1 The comoving number densities of SMBHs in different mass bins at redshift z — 6. The 
24 different models shown in the figure assume different parameter combinations as follows. The 
columns, from left to right, adopt / see( j = 10~ 3 , 10~ 2 , 10~', 1. The top row assumes a random 
binary spin orientation, and the bottom row assumes that BH spins are aligned with the binary's 
orbital angular momentum. In each panel, the time-averaged mass-accretion rates, in Eddington 
units, are assumed to be either /duty = 1 (black solid curves), /duty = 0.8 (blue dotted), and /duty = 
0.6 (green dash-dotted). The numbers in the upper-right corners represent the total mass density in 
SMBHs log 10 [p./(M;. ; Mpc~~ )] for each model. The red dashed line demarcates the abundance of 
z w 6 SMBHs with m > 1O 96 M already observed in the SDSS (adapted from ref. l250l ). 



with a Navarro-Frenk-White (NFW) profile H1771 . and a more cuspy baryonic com- 
ponent (with p oc r~ 2 - 2 , suggested by 3D simulations). At the time of a merger, the 
trajectories of kicked BH - ejections, or oscillations damped by dynamical friction 
- were followed explicitly by one-dimensional orbital calculations. 

Figure Q] shows that a robust conclusion for a model to produce enough z = 6 
SMBHs is that /d uty > 0.6 - namely the w 100 M stellar seed BHs must accrete 
near the Eddington rate nearly all the time. (Note that this value is excluded in 
the spherically symmetric case H1741 .) The initial BH occupation fraction also has 



10 Zoltan Haiman 

to be / S eed > 10~ 3 . Another interesting, and less intuitive conclusion, is that if the 
initial seeds are rare (/ see d = 10~ 3 — 10~ 2 ), then gravitational kicks do not have 
a big impact, and it makes little difference to the SMBH mass function whether 
spins are aligned or randomly oriented (this can be seen by comparing the bottom 
and top panels in Fig. [D. This is because the few "lucky" seeds that form earliest 
(at z > 25) have a chance to grow by > two orders of magnitude in mass before 
encountering their first merger. The masses of the two BHs at this first merger are 
then very unequal (q = M1/M2 < 0.01), making kick velocities too low to lead to 
ejection. It is important to emphasize, however, that the model trajectories for the 
kicked BHs assume spherical symmetry II 1561 |3T1 [80l 12501 . In gas-rich galaxies, 
most of the dynamical friction occurs due to the dense baryons at the center of the 
potential |250 82 1. In asymmetric potentials, the kicked BH does not return to the 
central region - its oscillations are damped less quickly fl8~Tl and the accretion rate 
onto the oscillating hole is also suppressed. 

An important additional issue is that in those models that satisfy the observational 
constraint on the SMBH abundance, the massive end of the SMBH mass function 
is extremely steep. This prediction is not surprising, as the most massive SMBHs 
reside in few x 10 12 M Q halos, which probe the 5<7 tail of the halo mass function at 
z = 6 (and there are indeed as 10 8 (!) times as many fewx 1O 9 M halos, which host 
~ 10 6 M e BHs). As a result, the total mass density in SMBHs with masses above 
> 10 5 M© BHs (shown by the numbers in the upper right corners in Figure Q]) are 
overpredicted by a factor of 10 2 — 10 3 . Note that these numbers indicate the mass in 
SMBHs that avoided ejection due to kicks, and remained in galactic centers (in some 
of the models, a significant fraction of the BHs are ejected and form intergalactic 
BHs; there is no obvious means to detect these 12501 ). The mass density of such 
nuclear SMBHs at z ~ can be inferred from the observed correlations between BH 
masses and host galaxy properties (such as the masses or velocity dispersions of the 
host halos; e.g. [64 1). The result is several x 10 5 M Q Mpc~ 3 ; furthermore, the expec- 
tation is that most (> 90%) of this mass was accreted well after z = 6 H218I . Some 
strong feedback is therefore needed to eliminate this significant overprediction. Pos- 
sible candidates for this are radiative feedback internal to halos, which maintains the 
"M — CJ relation" in ultra-high redshift, low-mass halos, or the termination of PopIII 
star formation, at redshifts as high as z ~ 20, due to Lyman Werner radiation ll89l or 
metal pollution l36l . 

Finally, it is worth emphasizing that the mass accretion rate corresponding to the 
Eddington limit - for the fiducial radiative efficiency of e = L/ihc 2 =0.1 for con- 
verting mass to radiation - would need to be exceeded only by a factor of a ~Tew to 
make the growth from stellar seeds much easier. Modestly exceeding the Eddington 
rate is theoretically certainly plausible (see below): density inhomogeneities can al- 
low radiation to leak out of low density regions while most of the accreting matter 
can be contained in high density regions. For example, magnetized radiation dom- 
inated accretion disks are subject to a "photon bubble" instability that nonlinearly 
appears to lead to strong density inhomogeneities (e.g. |[T8l ). Nevertheless, obser- 
vations have so far not revealed systems that sustain super-Eddington accretion for 
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extended periods; it would then still have to be explained why the z ~ 6 quasar BHs 
have this uniquely different behaviour. 



3.1.3 Accretion versus Mergers 

Mergers between halos can help build up the mass of individual black holes (with- 
out significantly changing the total mass of the population), provided that the central 
black holes in the halos coalesce rapidly. The mean accretion efficiency of ~ 10% 
inferred from comparing the local black hole mass density with the integrated quasar 
light suggests that accretion dominates at least the last e-folding of the black hole 
mass 1302||218l . Mergers may, however, be significant earlier on |9T|. In addition, 
uncertainties in the expected radiative efficiency of black hole accretion limit how 
accurately one can constrain the growth of black hole mass by mergers. For exam- 
ple, if the typical efficiency was sa 40%, as for a maximally rotating Kerr black hole, 
then the Eddington-limited mass accretion rate would be decreased correspondingly, 
and mergers could dominate black hole growth (on the other hand, note that multiple 
mergers would have a tendency to cancel the black hole spin; HI 121 ). In order for 
mergers to contribute significantly to the growth of individual black hole masses, 
stellar seeds must be present in large numbers, in the most of the numerous miniha- 
los that form at z > 15, down to small halo masses. 

The balance between growth through BH mergers and growth through gas accre- 
tion is indeed a key characteristic of any SMBH assembly scenario. For concrete- 
ness, consider possible merger histories for the z = 5.82 SDSS quasar SDSS 1044- 
0125 ([99|, the following arguments would be stronger for more luminous quasars 
at higher redshift). One can estimate the mass of the dark matter halo harboring the 
quasar by its abundance. SDSS searched a comoving volume of ~ 1 Gpc 3 to find 
each quasar. Assuming a duty cycle of a few times 10 7 years, one estimates that the 
dark matter halos corresponding to this space density have masses of 10 13 M Q (us- 
ing the halo mass function in ref. 11181 . the original Press-Schechter formula 111931 
would give a similar answer). A 10 13 M© halo at z = 6 typically has only ~ 10 
progenitors with circular velocities of v > 50 km s (the other progenitors being 
smaller). This implies that mergers can only help build up the black hole mass if 
seed black holes are present in progenitor halos with much smaller masses. A cutoff 
in the black hole mass function in halos with circular velocities below v = 50 km s _1 
would be justified if the cosmic ultraviolet background could suppress gas infall into 
smaller halos I56ll254lll78lll33l . However, one-dimensional gas collapse models 
with radiative feedback |f52]| have shown that this suppression is ineffective at red- 
shifts beyond z > 6. Thus, there is no obvious obstacle to forming seed black holes 
in halos down to v ~ 10 km s _1 (below this threshold, atomic H cooling becomes 
inefficient, and H2-photodissociation can be a limitation). 

In the illustrative models discussed in the previous section, which successfully 
reproduce the abundance of the x 10 9 M H SDSS quasar BHs, gas accretion accounts 
for the vast majority of the growth (in the sense that if the seed BHs were simply 
added together without any further accretion, the resulting total BH mass at z = 6 
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would be reduced by many orders of magnitude; see Table 3 in H250I ). However, 
in versions of these models in which BH growth is assumed to be self-regulating, 
accretion is much less important. Such models essentially describe the most heavily 
merger-driven scenarios possible, requiring accretion-driven growth of as little as a 
factor of a few. This is not surprising: placing a seed black hole in each arbitrarily 
low mass progenitor halo, with the same black hole mass to halo mass ratio as 
inferred for the SDSS quasars (Af./Afhaio ~ 10~ 4 ), could account for the observed 
black hole masses in quasars by z — 6, even without any gas accretion J9T1 . 

A further important unsolved question is whether halo mergers necessarily lead 
to black hole mergers at all (see, e.g., H1681 for a review). During a galaxy merger, 
the black holes sink via dynamical friction to the center of the galaxy and form a 
tight black hole binary in the nucleus. In normal galaxies with a stellar component 
in the nucleus, the black hole binary can continue to shrink by ejecting low-angular 
momentum stars that pass close to the binary (those in the "loss cone"). This process, 
however, clearly does not operate in the earliest stages of structure formation, when 
there are at most a few stars (if any) present in the merging mini-galaxies. Even 
at the later states, this process is inefficient, at least in spherical stellar systems, 
because the loss cone must be replenished by two-body relaxation. The black hole 
binary thus appears to stall and cannot coalesce even during a Hubble time [20! . 

Several ideas for circumventing this difficulty have been proposed. At later 
stages, in triaxial stellar populations, low-angular momentum orbits are populated 
much more efficiently because the stellar orbits can be chaotic; the resulting binary 
decay times are in many cases significantly less than a Hubble time, even if only a 
few percent of the stellar mass is on chaotic orbits (e.g., H301lll69lll94|[T58l ). In the 
earliest galaxies without large stellar populations, the coalescence of BHs must be 
facilitated by gas physics. If circumbinary gas is present and forms a thin accretion 
disk, then BH-disk interactions can drag the binary together, in a manner similar to 
Type II migration in planetary systems. The main difference from the planet case is 
that at least for nearly equal mass binaries, the secondary BH's mass will far exceed 
the mass of the disk, which slows down the migration (see |95l for a comprehensive 
discussion, and [76,8. 60| for examples of earlier work). A binary BH embedded in 
a spherical gas cloud is also facilitated by gaseous torques ll60l ; this case is much 
less well explored, but is likely to be more relevant to the earliest stages of the 
growth in the stellar-seed models, in halos without stars and with potential wells too 
shallow to support a thin disk. Finally, if SMBHs are brought together by successive 
halo mergers at a rate higher than the rate at which they can coalesce, then one or 
more of the BHs can be ejected out of the nucleus of the merger remnant by the 
slingshot mechanism 112091 . This could have implications for SMBH mass build-up 
in principle; in practice, more recent work on the dynamics of triple BHs indicate 
that ejections are relatively rare, and in the majority of cases, at least two of the BHs 
coalesce ifTTOllTl. 
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3.1.4 Super-Eddington Mass Accretion 

If mass is supplied to a black hole at m = M /MEdd 3> 1> the photons are trapped 
in the inflowing gas because the photon diffusion time out of the flow becomes 
longer than the time it takes the gas to accrete into the black hole (e.g. Ifl6l |2T1 ). 
The resulting accretion is thus not via the usual thin disk H2161 . but rather via a 
radiatively inefficient flow (RIAF). The luminosity is still set by the Eddington limit, 
but most of the gravitational binding energy released by the accretion process is not 
radiated away (being trapped in the flow). 

It is attractive to assume that the growth of SMBHs at high redshifts proceeds via 
such an optically thick, photon trapped accretion flow with m 3> 1 . Indeed, it would 
be a remarkable coincidence if the mass supply rate were precisely ~ MEdd (required 
for a thin accretion disk) during the entire growth of massive black holes. It is more 
likely that the mass supply rate is initially much larger in the dense environments of 
high redshift galaxies (ra>l) and then slowly decreases with time as the galaxy is 
assembled and the BH gains mass (e.g., J239l l40l ). 

Three-dimensional simulations for the cooling and collapse of gas into the 
first minihalos find that H2 cooling reduces gas temperatures to a fewxlOOK and 
produces a quasi-static contraction, with relatively low mass accretion rates of 
10~ 3 — 10~ 2 M yr~' (e.g. Q). This external mass-supply rate would still cor- 
respond to a super-Eddington growth rate for a BH with a mass of < 10 5 M { 7 ,. How- 
ever, there is another limitation: within a Kelvin-Helmholtz time of ~ 10 5 years, 
only a few x 100 M of material is accreted to the center (this is shown explicitly 
in Figure|Ui. Much more mass than this is then unlikely to be incorporated onto the 
central proto-star, before it settles to the main sequence. Radiative feedback from 
the proto-star (in the form of H2 dissociation, Lyman a radiation pressure, and ul- 
timately, photoionization heating) on the infalling envelope was found to limit the 
final mass of the star to » 140 M B249II165I . 

On the other hand, in more massive halos provided that H2 cooling can be dis- 
abled throughout the entire time of the collapse, the gas temperature is set by atomic 
H cooling and remains near 10 4 K. In a self-gravitating gas, the mass accretion rate 
is of order ~ c^/G <* r 3 ' 2 /G (e.g. 02341 ). Three-dimensional simulations have con- 
firmed this scaling (e.g. 1187112171 ). and in halos with r v ; r ~ 10 4 K, have found mass 
accretion rates of ~ 1M @ yr _1 112171 . As shown explicitly in Figure @] with this 
higher mass accretion rate, the mass that can be accumulated in the nucleus within 
a Kelvin-Helmholtz time is increased to 10 5 M . 

Theoretical models for the accretion on much smaller spatial scales (not resolved 
in the above simulations) imply that even if m ^> 1, only a small fraction of the mass 
supplied to the black hole actually reaches the horizon; most of it is driven away in 
an outflow (see, e.g., simulations of RIAFs B2481 12471 [l05l 1 1 1 41 1 1 961 : and analytic 
models [29, 30] 11971 ). The accretion rate onto a black hole thus probably cannot 
exceed ~ MEdd by a very large factor, even if the mass supply rate from larger radii 
is large (see 112161 for an early discussion of this point). 

The above discussion focuses on whether highly super-Eddington accretion is 
possible. The question of whether the Eddington limit for the luminosity can be 
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exceeded by a modest factor of ~ 10 is a bit more subtle. Magnetized radia- 
tion dominated accretion disks are subject to a "photon bubble" instability that 
nonlinearly appears to lead to strong density inhomogeneities (see, in particular, 
[9, 68, 17, 28, 18]). Density inhomogeneities allow super-Eddington fluxes from the 
accretion flow because radiation leaks out of the low density regions while most of 
the matter is contained in high density regions. Ref. |[T8l estimates that the Edding- 
ton limit can potentially be exceeded by a factor of ~ 10 — 100. This would allow 
much more rapid growth of black holes at high redshifts, circumventing the above 
arguments that seed black holes at z ~ 15 are required. Magneto-hydrodynamic 
(MHD) simulations of radiation dominated accretion flows have confirmed the rapid 
growth of unstable short-wavelength modes, with the development of large density 
variations. Inhomogeneities then allow the radiation to diffuse outward five times 
more rapidly than in a disk in hydrostatic equilibrium with no magnetic fields [260]. 
Explicit models for such slim, porous, accretion disks have been constructed re- 
cently l54l . In these models, when the external mass accretion rate is 10-20 xAife^, 
despite the presence of winds, a super-critical fraction, 2.6-3.8 MEdd> was indeed 
found to reach the central SMBH. 



3.2 Growth by Rapid Direct Collapse 

An appealing alternative idea is to produce, say, a 10 5 Mq SMBH "directly" - i.e. 
much faster than this would take under Eddington-limited accretion from a stellar 
seed. This would clearly be helpful to explain the high-redshift SMBHs. In this con- 
text, the crucial question is whether gas can accrete at a highly super-Eddington rate 
onto a black hole, i.e., with M^>M Edd , where Mgdd = ^L Edd /c 2 « 1.7MgM Q yr _1 
is the accretion rate that would produce an Eddington luminosity if accretion onto a 
black hole of mass 10 s M$ M Q proceeded with 10% radiative efficiency. If so, this 
could lead to rapid black hole growth at high redshifts. Constraints on BH seeds and 
their formation redshifts would therefore be much less stringent. 



3.2.1 Rapid Collapse of Gas in r vir > 10 4 K Halos 

SMBHs may form directly by the collapse of gas clouds at high redshifts, likely 
via an intermediate stage of a supermassive star or disk. The gas must not only 
shed angular momentum efficiently and collapse rapidly, but must also then avoid 
fragmentation. Whether fragmentation of the gas cloud into stars can be avoided is 
particularly questionable, in view of the large angular momentum barrier that must 
be overcome to reach small scales in a galactic nucleus (forming an SMBH through 
a dense stellar cluster is another option, discussed in the next section). 

The most promising locations for such rapid "direct collapse" are at the centers 
of halos with r v ; r ~ 10 4 K. In the past several years, many authors have sketched 
how gas may collapse rapidly, without fragmentation, in these halos. The essential 
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idea is that when contracting gas in a protogalactic nucleus becomes optically thick 
and radiation pressure supported, it becomes less susceptible to fragmentation and 
star formation. It is, however, unlikely that radiation pressure becomes important be- 
fore angular momentum does, implying that the gas forms a viscous accretion disk 
in the galactic nucleus (fragmentation before the disk forms can also be avoided 
if the forming fragments collide and "coalesce" before they can separate into dis- 
crete dense clumps 1127] ). On the other hand, if self-gravitating, the resulting disk is 
strongly gravitationally unstable and becomes prone to fragmentation and star for- 
mation (e.g., 1123311741 ). Whether this fragmentation can be avoided is unclear. One 
possibility is to stabilize the disk by keeping its temperature "warm" (i.e. T ~ 10 4 K, 
close to the virial temperature). This would fatten the disks (the scale height scales 
with the ratio of gas and virial temperatures); this scenario may be possible in a 
virtually metal-free, high redshift halo flT82] |37] |288] |203] . 

A suite of recent numerical simulations studied gas collapse in halos with r v ; r ~ 
10 4 K 12171 . It was found found that the gas in such halos, when collapsing in isola- 
tion, forms H2 efficiently, and (unfortunately) cools to temperatures of T ~ 300K. 
Although no fragmentation was seen, the cold gas (well below the virial temper- 
ature) is expected to ultimately fragment on smaller scales that have not yet been 
resolved [258 1. More importantly, even if fragmentation was successfully avoided, 
there is a problem: the cold gas was found to flow inward at relatively low veloci- 
ties, near the sound speed of ~ 2 — 3 km s , with a correspondingly low accretion 
rate of ~ 0.01 M yr _1 . 02171 speculate that this is explained by a series of weak 
shocks in the infalling gas, which prevent the gas from accelerating to large Mach 
numbers (this is similar to the behavior seen in three-dimensional simulations of the 
so-called "cold mode" of accretion in lower-redshift galaxies 11311 ). Ultimately, the 
slow infall velocities and cold temperatures produce conditions nearly identical to 
those in the cores of lower-mass minihalos (mentioned above); extensive ultra-high 
resolution simulations had concluded that the gas then forms a single ~ 100 M© star 
Il3l l35ll299l or perhaps fragments even further into several stars [258 2461 1771 fl 951 . 
rather than forming a supermassive star or BH. 

There have been at least three different ideas on how to avoid H2-cooling and 
to keep the gas warm. One is for the gas to "linger" for a sufficiently long time 
at 10 4 K that it collapses to a SMBH, even before H2 has a chance to reduce the 
temperature (H2 is kept dissociated by collisions before the temperature falls be- 
low ~ 4,000K). For a sufficiently high space- and column-density of neutral hy- 
drogen, the absorption of trapped Lyman a photons can be followed by collisional 
de-excitation, rather than the resonant scattering of the Lyman a photon, effectively 
trapping much of the cooling radiation. This could prevent the gas temperature from 
falling below ~ 8,000K, and lead to such lingering and to SMBH formation - anal- 
ogous to opacity-limited fragmentation in colder gas in the context of star formation 
12401 12 101 [1441 . 

Another possibility is that, even in the presence of significant cooling, angular 
momentum transport by gravitational instabilities, spiral waves, bars, etc., can drive 
a fraction of the gas to yet smaller scales in the galactic nucleus (e.g. f24l ). Ref. |58l 
argued that this was particularly likely to occur in rare low angular momentum dark 
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Fig. 2 Temperature evolution of a metal-free cloud, irradiated by a strong UV flux. The models 
solve for the chemical and thermal evolution, but assume a pre-imposed density evolution, based 
on the spherical collapse model. Various cases are shown, with UV intensities at the Lyman limit 
of 721 = 0, 1, 10, 100 and 10 3 , in the usual units of 10~ 21 erg cm~ 2 sr~' s -1 Hz 1 (solid and dashed 
curves; see the legend in the panel). Each blue dotted line corresponds to a different constant Jeans 
mass. The gas is heated adiabatically until a density of n « 10° — 10 2 cm~ 3 , at which H2-cooling 
becomes efficient and cools the gas to a few x 100 K. However, there exists a critical flux, with a 
value between /ji — 10 2 and 10 3 , above which H.2-cooling is disabled (adapted from ref. [ 186 |). 



matter halos because the disk could viscously evolve before star formation com- 
menced. A similar idea is that a small fraction of the gas, with low specific angular 
momentum, within the halo may collapse to the center without undergoing fragmen- 
tation IU38II1521 . It may help that even if most of the gas is initially converted into 
stars, stellar winds and supernovae will eject a significant amount of this gas back 
into the nucleus; some of this gas can eventually collapse to smaller scales 11221 . 

Finally, H2-cooling may be disabled if the gas is exposed to a sufficiently intense 
UV flux 7, either directly photo-dissociating H2 (in the Lyman-Werner bands near 
a photon energy of ~ 12eV) or photo-dissociating the intermediary H (at photon 
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energies > 0.76eV). Requiring the photo-dissociation timescale, f^iss x J > to be 
shorter than the ^-formation timescale, tf orm <x p~ 1 ! genetically yields a critical 
flux that increases linearly with density, J cnt oc p. In low-mass minihalos, the critical 
flux is low, 7 crit w 0.01 - 0.1 EOll [Bl EQl [T/TTFI . Since the gas in halos with 



?vir > 10 4 K can cool via atomic Lyman a radiation and loose pressure support, 
it inevitably collapses further. As a result, in these halos, the critical flux is high, 
/ cnt « 10 2 — 10 5 , depending on the assumed spectral shape (ref. 02171 ; see also refs. 
Ml 851 l37l who found similar, but somewhat higher values). The existence of this 
critical flux is illustrated in Figure [2] using a one-zone model, in which the density 
evolution is approximated by spherical collapse, and the gas is illuminated by a 
source with a black-body spectrum (with a temperature of T = 10 4 K, characteristic 
of a normal stellar population). Figure |3]shows the radial structure of a 1O 8 M halo, 
at the time of its collapse, when illuminated at various intensities, taken from three- 
dimensional simulations with the AMR code enzo. These profiles clearly show that 
when the UV flux exceeds a critical value, the core of the halo is prevented from 
cooling to low temperatures. 

The 3D simulations also provide an estimate of the mass of the central "object" 
(star or SMBH) that ultimately forms at the core of the halo, based on the mea- 
sured profile of the mass accretion rate. This is illustrated in Figure [4] In particular, 
when the flux exceeds the critical value, and the gas remains warm, the collapse is 
significantly delayed. However, when the gas ultimately does collapse, it accretes 
toward the center at the sound speed (c s ~ lOkm/s), leading to a mass accretion rate 
of M w 1M Q yr~ ' . This is much higher than in the case of cold (c s ~ 1 km/s) gas in 
halos with efficient H2 cooling (as mentioned above, the simulations find M °= c 3 , 
as expected in self-gravitating gas). 

Importantly, the critical flux is high - likely significantly exceeding the expected 
level of the cosmic UV background at high redshifts. Therefore, only a small subset 
of all r v ; r > 10 4 K halos, which have unusually close and bright neighbors, may see 
a sufficiently high flux. However, given the strong clustering of early halos, there is 
a sufficient number of these close halo pairs to account for the abundance of the rare 
z = 6 quasars |[5T) . A more significant challenge to this idea is that in order to avoid 
fragmentation, the gas in these halos must also remain essentially free of any metals 
and dust 11861 . This requirement could be difficult to reconcile with the presence of 
a nearby, luminous galaxies. 

Finally, an important point to emphasize is that the collapsing gas is optically 
thick, and the critical flux J a i t depends crucially on the details of self-shielding of 
the Lyman-Werner lines of H2. Since following radiative transfer in many dozens of 
lines is computationally expensive, existing works have employed various simpli- 
fying approximations. The simplest (and by far most commonly used) approach is 
to combine a simple power-law fitting formula, /shield = (A^/lO^cm -2 ) -3 ' 4 , for 
the H2 self-shielding factor 1331 with an estimate for an effective H2 column density 
A^h, (most often equated with the product of the local density and Jeans length). 



4 Here J denotes the specific intensity just below 13.6eV, in the usual /21 units of 
10~ 21 erg cm~ 
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Fig. 3 The results of adaptive mesh refinement (AMR) simulations of a primordial halo, with a 
total mass of a few xlO 7 M , collapsing at redshift z ~ 10—15, exposed to various UV background 
fluxes. The four panels show snapshots of the spherically averaged profile of the particle number 
density, gas temperature, e~ fraction and Hi fraction, at the time of the collapse of the core, for 
several different values of the UV background intensity J21 , as labeled. The existence of a critical 
flux, here with a value between J21 — 30 and 10 2 , above which tb-cooling is disabled, is evident 
(adapted from ref. 12171 ). 



These assumptions have recently been scrutinized in ref H289I , which solved radia- 
tive line transfer exactly, using post-processing of three-dimensional simulations. 
This showed, rather promisingly, that when self-shielding is treated more accurately, 
J cr ii is reduced by about an order of magnitude. Interestingly, this reduction comes 
from a product of three very different sources (each of which individually reduces 
the shielding by a factor of ~two): (1) a numerical inaccuracy of the power-law 
/shield formula, (2) the inapplicability of this fitting formula at the relevant, rela- 
tively high temperatures (> 10 3 K), where excited rotational levels of H2 are pop- 
ulated, and (3) the Jeans length yielding an overestimate of the effective average 
(over different sightlines) column density. The order-of-magnitude reduction in J 
is especially important, since the probability distribution of J, sampled by halos at 
z > 10, is very steep near J ~ 10 4 (see OP ). With the original high / cr j t value, it has 
been shown ||5T1 that only one in w 10 6 halos - only those with an unusually bright 
and close neighbour - will see a sufficiently high flux. The reduction of the 7 c rit,2i 
value will significantly increase the number of candidates for objects that can avoid 
H2-cooling and fragmentation, and makes this scenario much more viable. 
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Fig. 4 The local accretion time-scale ? acc as a function of the enclosed gas mass M eas , in the same 
halo depicted in Figure|3] illuminated with different intensities, as labeled. In the halos exposed to a 
supercritical flux (/21 = 10 2 and 10 3 ), the mass accretion rate, M « 1 M ( T :i yr _1 , is nearly 10 3 times 
higher than in halos whose gas cools via Ht {J21 < 10). At the center of the brightly illuminated 
halos, ss 10 M© of gas accumulates within a Kelvin-Helmholtz time of f» 10 5 M ( ;), possibly leading 
to the formation of a SMBH with a comparable mass (adapted from ref. 12171 ). 



3.2.2 The Ultimate Fate of the Gas 



Although the detailed evolutionary pathways are still not understood, a possible 
outcome of the above scenarios is the continued collapse of some gas to smaller 
scales in the galactic nucleus. As the gas flows in, it becomes optically thick, and 
the photon diffusion time eventually exceeds the inflow time. Radiation pressure 
dominates for sufficiently massive objects so that the adiabatic index is F « 4/3. 
Radiation pressure may temporarily balance gravity, forming a supermassive star or 
disk (SMS; e.g., Ml 11112771 ; see, e.g., M227I for a review and additional references 
to earlier work). The SMS will radiate at the Eddington limit and continue contract- 
ing. When the SMS is sufficiently compact (GM/Rc 2 w 10~ 4 M g for non-rotating 
stars), general relativistic corrections to the gravitational potential become impor- 
tant, and the star becomes dynamically unstable because its effective polytropic 
index is < 4/3. For masses < 10 M , thermonuclear reactions halt the collapse 
and generate an explosion (e.g., l66l ), but more massive objects appear to collapse 
directly to a SMBH (see M225B for a review; and, e.g., [231 206| for recent simula- 
tions). 

If the mass accretion rate is high (~ 1 M yr _1 ), the outer layers of the SMS 
do not have time to thermally relax, and a high-pressure core-envelope structure 
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may develop, dubbed a "quasi-star" |25l |23l [T9l . The envelope initially contains 
most of the mass, and the central BH embedded in the envelope can grow from 
this envelope; the key feature of this configuration is that the accretion is limited by 
the Eddington limit for the entire envelope, rather than just the BH. Interestingly, 
the mass accretion rate required for this model is comparable to that seen in three- 
dimensional simulations of T y „ > 10 4 K halos, provided again that the gas in these 
halos can avoid Hj_ cooling 112171 . 

Finally, a possibility for the gas is to ultimately fragment into stars, but not before 
it reaches very high densities. If the gas is metal-enriched, this scenario may be most 
likely, and will result in the formation of a dense and compact stellar cluster, which 
naturally evolves to form a SMBH [1861. We turn to this idea in the next section. 



3.3 The Formation of Black Holes in Stellar Clusters 

The negative heat capacity of self-gravitating stellar systems makes them vulnera- 
ble to gravitational collapse in which the core of the cluster collapses on a timescale 
t cc comparable to the two-body relaxation time of the cluster (27]. If core collapse 
proceeds unimpeded, the resulting high stellar densities can lead naturally to the 
runaway collisional growth of a single massive object which may evolve to form 
a black hole (as in the discussion of SMSs above). This process provides an addi- 
tional route for the direct formation of SMBHs at high redshifts (or, more likely, 
intermediate mass seeds). 

Early work suggested that the fate of stellar clusters depends sensitively on the 
number of stars in the cluster. 01471 and 111981 found that very dense massive star 
clusters (N > 10 6 — 10 7 stars) were required to have successful core collapse and 
runaway growth of a single massive object. In less massive clusters, core collapse 
was halted by binary heating, in which the cluster gains energy at the expense of bi- 
naries via three-body interactions H107II1131 . Successful core collapse also requires 
that t cc is shorter than the timescale for the most massive stars to evolve off the 
main sequence ( 01991 ; this requirement implies compact clusters < 1 pc in size). 
Otherwise, mass loss from evolved stars and supernovae prevents the core from 
collapsing (in much the same way as binary star systems can become unbound by 
supernovae). In principle, massive stars could evolve into stellar-mass BHs and form 
a dense cluster of stellar-mass BHs. In the context of high-redshift halo formation, 
this is naturally expected 01571 . given that the first stars in the first minihalos are 
believed to be massive. A dense cluster of stellar-mass BHs can, in principle, grow 
into a more massive IMBH by coalescence due to gravitational radiation; however, 
this process is effective only in large stellar systems that are found in present-day 
galactic nuclei II 1 841 . 

If the required number of stars was indeed as large as N > 10 6 — 10 7 , this would 
be bad news for early SMBH formation: in the cosmological hierarchy, > 10 8 M© 
halos - the smallest that could plausibly harbor such star clusters - are very rare 
(> 2.5<7 fluctuations) at z > 10. However, recent work has revived earlier ideas 
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that stellar clusters are subject to a "mass segregation instability" that makes even 
the relatively less massive clusters prone to forming black holes ([241 270, l22l ). 
Because massive stars in a cluster sink by dynamical friction toward the center (mass 
segregation), they invariably dominate the dynamics of the cluster core and can 
undergo core collapse on a timescale much shorter than that of the cluster as a whole 
(and on a timescale shorter than their main sequence lifetime). 01921 showed with 
N-body simulations that the resulting core collapse likely leads to runaway merger 
and formation of a single black hole, and |83l reached a similar conclusion for 
much larger N ~ 10 7 using Monte Carlo simulations. Including an explicit treatment 
of stellar collisions, the most recent Monte-Carlo simulations find that the central 
massive object that forms in core has a mass of ~ 10~ 3 of the whole cluster l75l . 

In the context of high-redshift BH formation, 11861 considered the cooling prop- 
erties of gas in r v ; r > 10 4 K halos. It is expected that the majority of such halos, 
when they are assembled, had already undergone some amount of star-formation. 
In this case, their gas would have at least a trace amount of metals. 11861 considered 
the case when such mildly polluted halos are exposed to a large UV flux, which 
dissociates H2, and initially prevents cooling. This allows the gas to contract to very 
high densities, without fragmenting initially. By following the thermal and chemical 
evolution of such low-metallicity gas, exposed to extremely strong UV radiation, 
Ml 861 found, however, that eventually, gas fragmentation is inevitable above a crit- 
ical metallicity, whose value is between Z cr » 3 x 10 _4 Z Q (in the absence of dust) 
and as low as Z cr » 5 x 1O~ 6 Z (with a dust-to-gas mass ratio of about O.OIZ/Z0). 
When the metallicity exceeds these critical values, an ultra-dense cluster (the den- 
sity at the time of fragmentation is n > 10 10 cm~ 3 ) of low-mass stars may form at 
the halo nucleus ( 11441 and 11491 argued for similar scenarios, in the central regions of 
a protogalactic disk). Relatively massive stars in such a cluster can then rapidly co- 
alesce into a single more massive object, which may produce an intermediate-mass 
BH remnant with a mass up to M < 10 3 M Q . 

The above processes provide a promising channel for the formation of IMBH 
seeds, which can grow via mergers and/or accretion to form SMBHs. For example, 
112721 and Ml 161 have incorporated such early black hole seeds into Monte Carlo 
simulations of the black hole merger histories. With reasonable prescriptions for the 
merging and accretion of black holes inside dark halos, these models can account 
for the observed evolution of the quasar luminosity functions at z < 5 and can serve 
for physically motivated extrapolations to high redshifts to describe the first AGN. 

It should be noted that there exist some observational evidence for IMBHs in the 
local universe. In particular, the presence of IMBHs with masses of order ~ 10 4 M Q 
are inferred from stellar kinematics in the globular clusters Gl ( 11711 ; note that the 
velocity dispersion profile itself does not require a BH |fl31 and the evidence comes 
from higher-order moments of the velocity distribution [f72l instead), (0 Cen Ml 801 
and in M15 ( 12631 , although this object can also be modeled without an IMBH; 
M262I ). Ultra-luminous X-ray sources in nearby galaxies (e.g., fl45l 11241 11721 l63l ) 
have also been interpreted as accreting IMBHs. While there are viable non-IMBH 
interpretations of these sources (e.g., 01321 H"8l ). the X-ray spectrum of one such 
source during the peak of an outburst implies the presence of a > 2000 M IMBH 
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(assuming that the luminosity is limited to 0.3 x L^a as in stellar-mass BH X-ray 
binaries in their hard state; 111231 ). 



3.4 Alternative Models 

Since both of the "standard" scenarios discussed above require some optimistic as- 
sumptions, it is interesting to consider some more exotic possibilities. 

It is commonly believed that the magnetic fields permeating galaxies such as 
the Milky Way arose by the amplification of a much weaker large-scale seed field. 
Weak primordial magnetic fields, with strengths of up to ~ InG, can be produced in 
phase transitions in the early universe, during inflation, or during the electroweak or 
QCD phase transitions. It has recently been shown that such a primordial magnetic 
field could produce a variant of the "direct collapse" scenario 12151 . In particular, if 
the field is tangled, then ambipolar diffusion will provide an efficient new mecha- 
nism to heat the gas as it collapses in protogalactic halos. If the field has a strength 
above B |> 3 (comoving) nG, the collapsing gas is kept warm (T ~ 10 4 K) until it 
reaches the critical density « cr it w 10 3 cm~ 3 at which the roto-vibrational states of 
H2 approach local thermodynamic equilibrium. This is illustrated explicitly by the 
thermal evolution of fluid elements shown in Figure [5] H2-cooling then remains in- 
efficient, and the gas temperature stays near ~ 10 4 K, even as it continues to collapse 
to higher densities. The critical magnetic field strength required to permanently sup- 
press H?-cooling is somewhat higher than upper limit of ~ 2nG from the cosmic 
microwave background (CMB). However, it can be realized in the rare > (2 — 3)(7 
regions of the spatially fluctuating B-field; these regions contain a sufficient number 
of halos to account for the z ~ 6 quasar BHs|j 

Another "exotic" idea is that the first PopIII stars may be powered by heat- 
ing by dark matter annihilation, rather than by nuclear fusion 12431 . Weakly in- 
teracting massive particles (WIMPs), can be such a heat source, as long as they 
reach sufficiently high density inside the first stars, and if the annihilation prod- 
ucts are trapped inside the star. Several authors have studied the impact of this 
additional heating mechanism on the structure and evolution of such "dark stars" 
EH] Hl5l |297] [2521 ESU EH] |65] EM). In particular, these stars can live much 
longer than "normal" PopIII stars, and because their radiation is soft, they can con- 
tinue to accrete gas, as long as the dark matter heating persists, and grow to masses 
of up to ~ 1O 5 M0. In fact, one of the challenges in these models is to explain why 
and how the growth of the star stops II261II651 . An interesting prediction is that these 
stars are bright, and should be detectable directly by JWST J65l . 



5 Because of the high magnetic Jeans mass, the magnetic pressure has significant dynamical effects, 
and can prevent gas collapse in halos with masses up to M > few x 10 10 M,t ; . These are ~100 times 
more massive than the DM halos in the "usual" direct collapse models discussed in § l3.2l above. 
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Fig. 5 The temperature evolution of a patch of the intergalactic medium is shown as it initially 
expands and then turns around and collapses to high density. The different curves correspond to 
different values of the assumed primordial magnetic field, as labeled. The gas evolves from the 
left to the right on this figure. The left panel shows the expanding phase, starting from an initial 
density of ss 100 cm~ 3 (corresponding to the mean density at redshift z ~ 800) and ending at 
the turnaround just below n = 10~ 2 cm~ 3 . The right panel follows the subsequent temperature 
evolution in the collapsing phase. The figure shows the existence of a critical magnetic field, with 
a value between 6 = 3 and 4 nG, above which Hi-cooling is disabled, and the gas temperature 
always remains near 10 4 K (adapted from ref. 12151 ). 



4 Observational Considerations 



In this section, we first review several recent observations and their implications for 
the formation of black holes at high redshifts. We then then speculate on how future 
observations may probe the assembly of high-z SMBHs. 
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4.1 Surveys for High Redshift Quasars 

The majority of the ~ 40 quasars known at z ~ 6 to date have been discovered in the 
SDSS. This is perhaps somewhat surprising, since the SDSS is a relatively shallow 
survey (with a magnitude limit of i ~ 22) capable of detecting only the rarest bright 
quasars at redshifts as high as z ~ 6. Nevertheless, the large solid angle searched for 
high redshift quasars (~ 8,000 square degrees) has yielded many such objects ll6D . 
Somewhat deeper surveys covering smaller areas (a few ^100 square degrees), 
such as the SDSS Deep Stripe [Qj]D, and the CFHQS (284) and UKIDSS Ifl46l 
surveys, have yielded many additional detections. The most important properties 
(for our purposes) of these sources are that they are probably powered by SMBHs 
as large as a few x 10 9 M and overall, they appear to be indistinguishable from 
bright quasars at moderate (z ~ 2 — 3) redshifts, with similar spectra and inferred 
metallicities (e.g. 11211 . In addition, a large reservoir of molecular gas is already 
present, even in the most distant sources 1 2781 12791 12801 . 

Despite the overall similarities, there are some tentative distinctions between 
these z > 6 quasars and their z ~ 2 — 3 counterparts. First, there is evidence for 
increasing Eddington ratios toward high redshift (see below). Second, quasar clus- 
tering has been found to strongly increase with redshift H229I . The observed clus- 
tering strength can be used to infer to quasar lifetimes Il93l 11611 , and implies that 
the duty cycle of bright quasar activity increases significantly toward high redshifts, 
to near unity by z ~ 6 [220|. Finally, two of the ~ 6 quasars have no detectable 
emission from hot dust [120|. There are no known examples for such apparently 
hot-dust-free quasars at low redshift; this result therefore suggests that at least these 
two quasar BHs may have been caught at a young age of the evolution of their host 
galaxies (i.e., there was insufficient time for a dusty torus to form in the nucleus). 

With these exceptions, the high-z SMBHs and their surroundings appear as "fully 
developed" as their lower redshift counterparts, despite the young age (< 10 9 years) 
of the universe at z > 6. These rare quasars are likely harbored by massive (~ 
10 12-13 M Q ) dark matter halos that form out of 4 — 5(7 peaks of the fluctuating 
primordial density field. The large halo mass follows directly from the space den- 
sity of these sources ( |99l ; another method to confirm the large halo masses is to 
study the expected Lyman a absorption signatures of cosmological gas infall onto 
such massive halos, as proposed by ll37l ). Indeed, the environment and dynamical 
history of an individual massive dark matter halo at z ~ 6 and z ~ 3 can be simi- 
lar; it is their abundance that evolves strongly with cosmic epoch. This is broadly 
consistent with the observations: the bright z ~ 6 quasars look similar to their z ~ 3 
counterparts, but their abundance is much reduced (by a factor of ~ 40). 

The fact that these quasars are so rare has important implications. First, they are 
likely to be the "tip of the iceberg" and accompanied by much more numerous popu- 
lations of fainter quasars at z > 6. The slope of the luminosity function is expected to 
be very steep at i ~ 22, and so pushing the magnitude limits further in future surveys 
should prove rewarding. The most direct constraints on this slope are from combin- 
ing the CFHQS sample with the more luminous SDSS sample (yielding a total of 
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40 quasars between redshifts 5 .74 < z < 6.42; 12841 ), and from gravitational lensing 
ll46l [292 293 204|. Combining source counts and lensing yields the strongest limit 
of — dlog<P/dlogL < 3 J291 II . Second, the steep slope of the dark halo mass func- 
tion implies that the masses of the host halos can be "measured" from the abundance 
quite accurately (see discussion in § I4.21 i. Conversely, since small changes in the as- 
sumed host halo mass results in large changes in the predicted abundance, large 
uncertainties will remain in other model parameters. In this sense, fainter, but more 
numerous quasars (or lack thereof) can have more constraining power for models 
that relate quasars to dark halos. 

The most striking feature of the SDSS quasars, however, is the large black hole 
mass already present at z ~ 6. In the rest of this section, we critically assess whether 
the inferred large black hole masses are robust. 

The masses of the black holes powering the SDSS quasars are inferred by as- 
suming that (1) they shine at the Eddington luminosity with a bolometric correction 
identical to that of lower redshift quasars (this is justified by their similar spectra), 
and (2) they are neither beamed nor gravitationally lensed (both of these effects 
would make the quasars appear brighter and allow lower BH masses). These as- 
sumptions lead to black hole masses M. « (2 — 6) x 10 9 M e for the z > 6 quasars 
known to date. These are reasonable assumptions, which have some empirical jus- 
tification. 

The hypothesis that the quasars are strongly beamed can be ruled out based on 
their line/continuum ratio. If the quasar's emission was beamed into a solid angle 
covering a fraction, /, of 4n, it would only excite emission lines within this cone, re- 
ducing the apparent line/continuum ratio by a factor /. However, the SDSS quasars 
have strong lines. ||90l found that the line/continuum ratio of the z = 6.28 quasar 
SDSS 1030+0524 is about twice that of the median value in the SDSS sample at 
z > 2.25 [264j. This argument, applied to the Mg II line of the z = 6.41 quasar 
SDSS Jl 148+5251 1128511 yields a similar conclusion. 

Another important uncertainty regarding the inferred black hole masses is whether 
the SDSS quasars may be strongly magnified by gravitational lensing. The optical 
depth to strong lensing along a random line of sight to z ~ 6 is small (~ 10~ 3 ; e.g., 
Ill34|[l4l ). Nevertheless, magnification bias can significantly boost the probability 
of strong lensing. If the intrinsic (unlensed) luminosity function at z ~ 6 is steep 
and/or extends to faint magnitudes, then the probability of strong lensing for the 
SDSS quasars could be of order unity [46 292 293 1). The overwhelming majority 
(more than 90%) of strong lensing events would be expected to show up as multi- 
ple images with separations at least as large as 0.3". It is difficult to produce strong 
magnification without such multiple images, even in non-standard lensing models 
(allowing ellipticity and/or external shear H130I ). However, deep optical observa- 
tions (e.g. with the Hubble Space Telescope) of the highest redshift quasars show 
no signs of multiple images for any of the z > 6 sources down to an angle of 0.3" 
12041 12861 fl 



6 The highest-redshift known lensed quasar is at z = 4.8 and was discovered serendipitously in 
SDSS, initially flagged as a galaxy due to the strong blending of one of the quasar images with a 



bright galaxy 11641 . 
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Finally, whether or not the SDSS quasars are shining at the Eddington limit is 
difficult to decide empirically. Several authors J285] J267] ED EQl |268] |269j IT391 



283| have estimated Eddington ratios in samples of high redshift quasars, using 
observed correlations between the size of the broad line region and the luminosity 
of the quasar (the correlation is calibrated using reverberation mapping of lower 
redshift objects; e.g. 1 128 2661 11361 ). Values range from sa 0.1 to > 1 ; in particular, 
L/LEdd is typically found to increase with redshift, and approaches unity for the 
z > 6 quasars. 

Inferences about Eddington ratios at high redshifts can also be made by utiliz- 
ing models of the quasar population as a whole. Such models typically assume the 
Eddington luminosity at higher redshifts, where fuel is thought to be readily avail- 
able [239 85 1 . Numerous semi-analytic models for the quasar population (see, e.g., 
ref. 112181 and references therein) have found that Eddington ratios of order unity 
during most of the growth of the black hole mass also yield a total remnant SMBH 
space density at z = that is consistent with observations. Ref. (43], and, more re- 
cently, 01731 have self-consistently modeled accretion and radiative feedback onto 
an individual quasar BH, and found that (provided fuel is available) the luminosity 
is near the Eddington value during the phases when the quasar is on. Despite these 
arguments, one cannot directly rule out the possibility that the SDSS quasars shine 
at super-Eddington luminosities (theoretically, this is possible, as discussed above). 
We emphasize that if this were true and the masses were lower than 10 9 M , then 
the SDSS quasars would have to be luminous for only a short time: maintaining the 
observed luminosities for > 10 7 years with a radiative efficiency of e = L/rhc 1 = 0.1 
would bring the black hole masses up to values of 10 9 M anyway. 



4.2 Local Black Holes as Fossils 

As mentioned above, SMBHs appear ubiquitous in local galaxies, with their masses 
correlating with the global properties of their host spheroids. Several groups have 
noted the broad natural implication that the formation of the SMBHs and their host 
spheroids must be tightly linked (see, e.g., 12181 ). Various independent lines of ev- 
idence suggest that spheroids are assembled at high redshifts (z ~ 2; see ||39l for 
the age determinations from the Sloan sample and references to older work), which 
would be consistent with most of the SMBH mass being accreted around this red- 
shift (coinciding with the peak of the activity of luminous quasars). Indeed, starting 
from the age distribution of local early-type Galaxies, one can reconstruct the cos- 
mic evolution of the quasar luminosity function to within observational errors, using 
the most naive set of assumptions (namely that the formation of stars and the assem- 
bly of the nuclear SMBHs track one another, with the SMBH radiating at a constant 
/Edd ~ 0.3, and that the M, — a correlation does not evolve with redshift [94, 12191 ). 
This then has the unwelcome (but unsurprising) implication that the local SMBHs 
may contain little direct evidence of the formation of their seeds at z > 6. Indeed, 
it seems most plausible that the observed tight correlations, such as between M. 
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and a, are established by a feedback process which operates when most of the 
black hole mass is assembled. However, an upside of this hypothesis is that — 
with the identification of a specific feedback mechanism — physically motivated ex- 
trapolations can be made toward high redshifts. Also, while relative massive local 
SMBHs (~ 10 9 M Q ) have undergone many mergers, those with the lowest masses 
(~ 10 6 M Q ) are more likely to have avoided mergers. Therefore, the low-mass end 
of the M. — a relation could be a probe of high-z SMBH formation models. To be 
more explicit: if only a small fraction of high-z halos are seeded with BHs, mas- 
sive galaxies will have undergone many mergers, and will have a nuclear SMBH at 
present (i.e. / see d ~ 1)- On the other hand, many low-mass galaxies may still have 
no BHs and thus the dwarf-galaxy population can have / see d <S 1 116711 . Likewise, 
the direct-collapse models produce SMBHs whose masses are initially well above 
the M. — a relation, which can be a diagnostic of such models [2501. Such an 'up- 
ward curvature in the low-mass end of the relation could indeed be preserved all the 
way down to z = 02741 . 

More generally, whether the local M. — a relation holds at higher redshifts, both 
in normalization and in slope (as discussed by several authors), and also in range 
(which has received less attention, but see 111791 ). are interesting observational ques- 
tions. The highest redshift SDSS quasars do appear to satisfy the M, — a relation of 
the local SMBHs, at least approximately. If M. is estimated assuming the Eddington 
luminosity, and a from the CO line-width, then, at a given a, high-z quasars have 
BH masses a factor of ~ 4 larger than local galaxies, although this could still be 
partly caused by selection (with high-z quasars preferentially viewed face-on) HI 091 
and by the observations probing a smaller (inner) fraction of the DM halo. We also 
note that if a is estimated from the circular velocity of the host dark matter ha- 
los with the right space density (e.g., J99l ), then the SDSS quasars are within the 
scatter of the M. — a relations of ll70l and also of J64). As explained in § 14.11 the 
halo mass inference is reasonable. The determination of the halo mass and circular 
velocity from the observed abundance of quasars is also more robust than it may 
at first appear. This is because, despite the dependence on the poorly known duty 
cycle, the halo mass function is exponentially steep for the massive M ~ 10 13 M 
halos at z ~ 6; therefore, the dependence of the inferred halo mass on the duty cycle 
(and other uncertainties in the estimated halo abundance) is only logarithmic. The 
weakest link in the argument is associating the spheroid velocity dispersion with 
the circular velocity of the dark matter halo. There is evidence [64] of a correlation 
between M. and o\ with the velocity dispersion measured in the dark matter dom- 
inated region of SMBH host galaxies; this establishes a direct link to the dark halo 
and puts the above argument on somewhat firmer ground (although there are still 
large errors in the inferred correlation, depending on the halo profile one adopts to 
convert the measured circular velocity to total halo mass). 

The (tentative) evidence that high-redshift AGN do not strongly deviate from 
the local M. — a relation further supports the idea that the formation of SMBHs 
and their host galaxies must be tightly coupled by cosmology-independent physi- 
cal processes (since the SDSS quasars are the rare peaks that have already formed 
at z ~ 6 instead of at z ~ 2). Besides the slope and normalization of the M, — a 
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relation, the range (of masses and velocity dispersions) over which observed galax- 
ies satisfy this relation has to match up between low and high redshifts. In par- 
ticular, the largest black holes observed at high redshifts have inferred masses ap- 
proaching M. ~ 10 10 Mq, These should also exist at low redshifts, but have not 
yet been discovered. In the SDSS, the galaxy with the velocity dispersion record 
has <r = 444 km s -1 12081 ). whereas a naive application of the local M. — a re- 
lation would predict the presence of <y > 700 km s -1 galaxies 111791 . This puzzle 
is alleviated somewhat by the scatter in the relation; it is likely fully resolved by 
the realization that the M. — a relation has 'curvature', with BHs in the largest 
galaxies preferentially more massive than the power-law M, — a relation would 
predict 01451 . 

There have been several suggestions in the literature for the nature of the dynami- 
cal coupling between the formation of the black hole and its spheroid host. The most 
promising is radiative or mechanical feedback from the SMBH on the gas supply in 
the bulge. The essential idea (going back to [236 1) is that when the black hole in the 
center of the galaxy grows too large, its outflows and radiation unbind the gas in the 
bulge or in the disk, quenching further black hole growth via accretion and further 
star formation. Competition with star formation for the gas supply may also play a 
role |50lll49| . Note that these mechanisms can readily work at any redshift. 

There are several alternative possibilities for the origin of the M. — a relation, 
which include: (1) filling the dark matter loss cone 111 881 . In this model, the growth 
of the SMBH occurs first through the accretion of collisional dark matter particles, 
and subsequently through the scattering of these particles into orbits that are then 
perturbed to pass sufficiently close to the black hole's Schwarzschild radius to be 
captured. This model runs into difficulties with the so-called Soltan argument; since 
the SMBHs are fed mostly dark matter rather than gas, there is no associated radia- 
tion. (2) Direct capture of stars on high eccentricity orbits by the SMBH |306l ll69ll . 
This model has a similar problem because black holes more massive than > 10 8 M 
do not tidally disrupt stars, so there is again no radiative output associated with the 
black hole growth. (3) Stellar captures by the accretion disk feeding the hole [ 175 1. 

Solving the puzzle of the origin of the M. — a relation will have important impli- 
cations for high-z SMBHs: in particular, it will generally determine how the relation 
evolves with redshift. 



4.3 The Future 

In this section, we briefly summarize the possibility of probing the continuum and 
line emission from AGN beyond the current redshift horizon of z ~ 6. This discus- 
sion is necessarily based on models for how the BH population evolves at z > 6. 
Such models can be constructed by assuming that SMBHs populate dark matter 
halos, e.g., in accordance with the locally measured M. — a relation (or an extrap- 
olation of the relation to higher redshifts). The relation appears to hold, at least to 
within a factor of a few, for z ~ 3 quasars (this is based on using the H/3/OIII lines 
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as proxies for black hole mass and a, respectively; e.g., [232|), and also at z ~ 6 
(see ij l4.ll ). No doubt the observational constraints will improve as both black hole 
masses and velocity dispersions are measured in larger samples of distant quasars. 
Correspondingly, extrapolations to high redshifts will be more reliable as the feed- 
back processes that regulate black hole growth are better understood. Here we sum- 
marize predictions from the simplest models. 



4.3.1 Broadband Detections 

Predictions for the number counts of high redshift AGN have been made using sim- 
ple semi-analytic models for the near-infrared 11961 and in the soft X-rays ll98l . In 
these early models, the quasar black hole was assumed to have a fixed fraction 
~ 10~ 4 of the halo mass, shine at the Eddington luminosity, and have a duty cy- 
cle of bright activity of t q ~ 10 6 years. 

In such models, the surface density of sources is very high in the optical/near- 
infrared bands, even at z ~ 10. For example, in the 1 — 5/im band, the ~ InJy sen- 
sitivity of the James Webb Space Telescope (JWST) will allow the detection of an 
~ 10 5 Mo black hole at z — 10 (provided that the black hole shines at the Eddington 
limit with a standard template spectrum ||59l ). Surface densities as high as several 
sources per square arcminute are predicted at this threshold from z > 5, with most of 
these sources at z > 10 11971 . We note, however, that these predictions are very sen- 
sitive to the assumed duty cycle of bright activity. For example, for tq ~ 10 7 years, 

5/3 

or M, <* Af h ^ , the z ~ 10 counts can be smaller by a factor of 10-100 (depend- 
ing on what redshift-dependence is assumed for the above scaling relation between 
black hole and halo mass at high redshift; see |96 84 295] for related discussion). 
It would also be interesting to detect the host galaxies of ultrahigh redshift AGN, 
which should be feasible with JWSTs sensitivity. If the galaxies occupy a fair frac- 
tion (~ 5%) of the virial radius of their host halos, then a large fraction (> 50%) of 
them can potentially be resolved with JWSTs planned angular resolution of ~ 0.06" 
ll97l [141 . The Large Synoptic Survey Telescope (LSSlO), with a planned capability 
of going ~ 5 magnitudes deeper than SDSS in a ~ 3 times larger solid angle, will be 
an ideal instrument for studying high-redshift quasars in the optical/near-infrared. 

In the soft X-rays, the 0.5 — 2 keV flux of 2.5 x 10~ 17 ergs cm~ 2 s reached 
in a 2 Ms exposure of CDF-North corresponds to a larger (~2x 10 7 M Q ; see 
Figure 1 in 11981 ) black hole at z = 10, but nevertheless, thousands of sources are 
predicted at z > 5 per square degree, and tens per square degree at z > 10. This 
would imply that tens of z > 5 sources should have been detectable already in the 
CDFs, whereas only a handful of potential candidates, and no confirmed sources, 
have been found. In revised models with longer quasar lifetimes and thus a steeper 
scaling of M. with Mhaio, these numbers can be sharply decreased ||97] [84], which 
can bring the expected counts into agreement with current non-detections H295I . 
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The radio sensitivity of the extended Very Large Array and other forthcom- 
ing instruments (e.g., Allen Telescope Array and Square Kilometer Array) is also 
promising for detecting AGN beyond z ~ 6. Using the updated scaling of black hole 
mass with halo mass and redshift [295] and assuming the same radio-loud fraction 
(~ 10%) as at lower redshifts, a simple model predicts that ^ten 10/xJy sources per 
square degree should be detectable at 1 — 10 GHz 11011 . The identification of these 
quasars is a challenge, but should, in principle, be feasible with deep optical/IR ob- 
servations. To date, only two such radio-selected quasars have been identified (in 
optical follow-ups of sources in the FIRST and VLA radio catalogs; 11631 13041 ). 
which falls below the expectations from the simple model by a factor of several. 

In addition to direct detection of AGN at very high redshifts, it may also be 
possible to detect lower mass seed black holes at comparable redshifts (or higher). 
In particular, a plausible model for gamma-ray bursts (GRBs) invokes accretion onto 
a newly formed ~ 10 M black hole (the collapsar model; e.g., 1 290 1). Swift has now 
detected four bursts beyond z > 6: GRB090429B at z = 9.4 gT), GRB090423 at z = 
8.2 EUIHD, GRB080913 at z = 6.7 (79], and GRB050904 at z = 6.3 |[T29lfT04l . 
for which the afterglow emission has also been measurable; such afterglows should 
remain detectable in the infrared out to z ~ 20 11411 l42l . Their detection and the 
characterization of their spectrum and light-curve would open up a new probe of 
black hole formation and evolution at high redshifts. 

In summary, model predictions for the continuum emission of z > 6 AGN are 
very sensitive to how one extrapolates the M, — Mhaio relation to z > 6. However, 
this should be viewed as "good news": (1) large numbers of detectable AGN at these 
redshifts are certainly possible, and (2) their detection will put strong constraints on 
models for the origin and evolution of the black hole population. 



4.3.2 Emission Line Measurements 

The strongest recombination lines of H and He from 5 < z < 20 AGN will fall 
in the near-infrared bands of JWST and could be bright enough to be detectable. 
Specific predictions have been made for the source counts in the Ha emission line 
11811 and for the three strongest Hell lines 1183112571 . The key assumption is that 
most of the ionizing radiation produced by the miniquasars is processed into such 
recombination lines (rather than escaping into the IGM). Under this assumption, the 
lines are detectable for a fiducial 10 5 M Q miniquasar at z — 10. The Lya line is more 
susceptible to absorption by neutral hydrogen in the IGM near the source but should 
be detectable for bright sources that are surrounded by a large enough HII region so 
that Lya photons shift out of resonance before hitting the neutral IGM |4~D . If the 
Lyman a emission is scattered off expanding shells of material (as expected from 
galactic winds), this will further shift the photons away from resonance, and make 
the emission line more detectable 1531 . We also note that in the "trapped" Lyman 
a model for direct collapse mentioned in § 13.21 the Lyman a emission ultimately 
emerging from the collapsing halo (before forming the SMBH) appears detectable 
with JWST, as well - as a low-surface brightness diffuse blob 111431 . 
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The simultaneous detection of H and He lines would be especially significant. 
As already argued above, the hardness of the ionizing continuum from the first 
sources of ultraviolet radiation plays a crucial role in the reionization of the IGM. 
It would therefore be very interesting to directly measure the ionizing continuum 
of any z > 6 source. While this may be feasible at X-ray energies for exceptionally 
bright sources, the absorption by neutral gas within the source and in the interven- 
ing IGM will render the ionizing continuum of high redshift sources inaccessible to 
direct observation out to lfim. This is a problem if the ionizing sources are black 
holes with M < 10 8 M at z ~ 10 (easily detectable at wavelengths red-ward of 
redshifted Lya in the near-infrared by JWST, but too faint to see in X-rays). The 
comparison of Ha and Hell line strengths can be used to infer the ratio of Hell to 
HI ionizing photons, Q = iV-™ /N™ ■ A measurement of this ratio would shed light 
on the nature of the first luminous sources, and, in particular, it could reveal if the 
source has a soft (stellar) or hard (AGN-like) spectrum. Note that this technique 
has already been successfully applied to constrain the spectra of sources in several 
nearby extragalactic HII regions J69l . Lyman break galaxies at z ~ 3 also appear 
to have unusually strong (for a normal stellar population) He 1640 emission line; 
however, the lack of X-rays rule out an AGN explanation (the observation could be 
explained instead by the presence of PopIII stars in these galaxies 11221 or by an 
unusual abundant population of Wolf-Rayet stars 0321 that can produce the He 1640 
line while avoiding an overproduction of metal lines). 

Provided the gas in the high redshift AGN is enriched to near-solar levels, several 
molecular lines may be visible. In fact, CO has already been detected in the hosts of 
the most distant quasars H278I 12791 l280l . The detectability of CO for high redshift 
sources in general has been considered in simple theoretical models 12371 . If AGN 
activity is accompanied by a star formation rate of > 30 M / yr, the CO lines are 
detectable at all redshifts z = 5 — 30 by the Millimeter Array (the redshift indepen- 
dent sensitivity is due to the increasing CMB temperature with redshift), while the 
Atacama Large Millimeter Array (ALMA) could reveal even fainter CO emission, 
and other C and O lines in emission, providing spatially resolved images [211 1. The 
detection of these molecular lines will provide valuable information on the stellar 
content and gas kinematics near the AGN. 



4.3.3 Gravitational Waves 

The most direct observational constraints on the SMBH assembly at z > 6, with 
especially clear distinctions between the "stellar-mass seed" and "direct-collapse" 
scenarios, may come from detecting the gravitational waves produced during the 
SMBH mergers. The Laser Interferometer Space Antenna (LISA) is expected to be 
able to detect mergers of SMBHs in the mass range ~ (1O 4 -1O 7 )M /(1 +z) with 
high S/N out to z ~ 30 |i2|. Binary spins and BH masses is expected to be mea- 
sured with high precision up to z ~ 10 |265|, especially if spin precession II 1421 and 
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higher-order harmonics of the waveform [11661 are included in the analysis^! Many 
authors have computed the expected LISA event rate from high-redshift SMBH 
merger population models in a range of plausible models. The published estimates 
(OSSISSllIHIIISIinillHET]; see a review in flOl ). even at lower redshifts, 
vary by orders of magnitude, from ~ 1 to as high as ~ 10 4 yr _1 ; there is a large 
range even among models that are explicitly calibrated to fit the evolution of the 
quasar luminosity function [11511 . 

Closest to the present context of the growth of SMBHs at z > 6 are the Monte- 
Carlo merger tree models in ref. [J250j| (discussed in § l3.1.2l above). These models 
coupled the merger trees with the orbits of oscillating kicked BHs, to predict detec- 
tion rates for LISA. They have surveyed a wide range of candidate assembly models, 
including those with rare, massive seeds, or through ultra-early production of nu- 
merous Pop-Ill remnant seeds. As mentioned above, in the latter model, seed BHs 
need to stop forming below a redshift z cut ~ 20, in order to avoid overproducing 
10 6 M BHs. 

The simplest SMBH assembly scenarios, which have constant accretion rates, but 
in which BH seed formation stops abruptly at some redshift, and which meet con- 
straints at both the high-mass and low-mass end of the z = 6 SMBH mass function, 
predict negligibly low LISA event rates. The reason for this pessimistic conclusion 
is as follows: in these models, the BHs that grow into the most massive, highest- 
redshift quasar-SMBHs accrete at the same (exponential) rate as all the other BHs, 
typically resulting in a vast overproduction of massive (in ~ 10 6 M©) holes. In order 
to offset this overproduction, seeds must be made very rare, and this diminishes the 
LISA rates. It is difficult to envision a scenario for high (> 10 per year per unit red- 
shift) detection rates unless a vast number of SMBHs in the 10 5 ~ 7 M© range lurk in 
the universe at all redshifts, which the current electromagnetic surveys have missed. 

A different class of models, which successfully build the z ~ 6 quasar BHs, are 
those in which the SMBH masses are self-regulated by internal feedback, to always 
maintain the M, — a relation. These models can evade this constraint, and produce 
LISA rates as high as 30 yr _1 . The key difference in these models with higher LISA 
rates is that the SMBH growth is driven by a large number of seed BHs and far 
lower gas accretion rates than those required in the constant-accretion models. The 
majority of the LISA events occur at z ~ 6 and in the low end (10 3 — 10 4 M ) of 
LISA'S mass range for detection. 

Also, for these models, the ejected BH mass density can exceed that of the galac- 
tic BH population at z = 6. Most ejected holes are expected to have low masses (still 
similar to the original seed mass), but an ejected BH can be as massive as ~ 10 8 M Q 
if large recoil velocities are allowed (e.g. if spins are not always aligned with the 
orbital angular momentum of the binary). 



8 As this paper was being written, NASA announced a decision to withdraw from the LISA exper- 
iment. The European Space Agency is continuing to consider a redesigned version of LISA, with 
a smaller budget, and a launch date of approximately 202 1 -2022. Given the very high S/N ratios 
forecast for the original version of LISA, the redesigned "LISA-lite" mission should still be able 
to detect low-mass SMBHs out to high redshifts. 
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Using similar "merger tree" massive black hole formation models, 02141 ana- 
lyzed the predicted mass- and redshift-distribution of LISA events. These models 
have input assumptions similar to the "M. — c" models in 1125011 . but with varying 
initial seed masses. These predict a handful of detectable events at z > 10. The raw 
total event rates in the two models are very similar. However, the mass-distribution 
of events is different (low-mass mergers are missing in the 'heavy seed' model). 
Another key diagnostic between the 'heavy' and 'stellar-mass' seed models is the 
mass ratio of the BHs in these detectable events: while the former models predict 
near-equal mass mergers, in the latter case, one of the merger partners typically have 
time to grow, resulting in typical mass ratios of q = 0.1 — 0.2. 

It is worthwhile to note that essentially all of the work on the gravity wave signal 
from black hole-black hole in-spiral has assumed efficient (nearly instantaneous) 
mergers. Stellar-scattering and gas can help drive BHs together on large scales, 
which can affect detections of the most massive nearby SMBHs by Pulsar Tim- 
ing Arrays 111351 and extreme mass ratio inspiral events by LISA 13031 . However, 
SMBH-SMBH coalescences, when they are in LISA'S frequency window, are well 
within the rapidly merging purely gravitational wave-driven regime, at least if the 
circumbinary gas forms a thin disk 



5 Conclusions 

In this review, we have summarized theoretical ideas and observational constraints 
on how massive black holes form at the centers of the earliest protogalaxies, and 
how such black holes grow via accretion and mergers to give rise to the observed 
population of black holes at z > 6 and in the local and moderate redshift universe. As 
this review shows, this remains a poorly understood but rich and important problem. 
Perhaps the most direct way of probing the role of mergers in black hole assembly 
and evolution at z > 10 is via their gravity wave signatures, which will hopefully be 
feasible with the redesigned version of LISA being considered by ES A. 

In addition to being of intrinsic interest for understanding the AGN phenomena, 
sources of gravity waves, etc., there is strong evidence that the formation and evo- 
lution of black holes is coupled to the formation and evolution of the host galaxy in 
which the black hole resides (e.g., the M. — a relation), and thus to the cosmological 
formation of nonlinear dark matter structures (i.e., the dark halos surrounding these 
galaxies). We anticipate that this will remain a growth area of research in the com- 
ing years, with continued rapid progress on both the observational and theoretical 
fronts. 
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